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Abstract 

Braneworld cosmology has several attractive and distinctive features. For instance the effective 



fN| . equation of state in braneworld models can be both quintessence- like (wq > —1) as well as phantom- 

>. 

like (wq < —1). Models with wq > —1 {wq < —1) are referred to as Brane2 (Branel) and 

o 

correspond to complementary embeddings of the brane in the bulk. (The equation of state in 



Branel can successfully cross the 'phantom divide' at w = — 1.) In this paper we compare the 
predictions of braneworld models to two recently released supernova data sets: the 'Gold' data 



Q-f (Riess et al, 2004) and the data from the Supernova Legacy Survey (SNLS) (Astier et al, 2005). 

We also incorporate the recent discovery of the baryon acoustic peak in the Sloan Digital Sky 
Survey (Eisenstein et al., 2005) into our analysis. Our main results are that braneworld models 
satisfy both sets of SNe data. Branel (with wo < —1) shows very good agreement with data for 

X' 

values of the matter density bounded from below: f^om ^ 0.25 (Gold) and f^om ^ 0.2 (SNLS). On 



the other hand Brane 2 (with wq > —1) shows excellent agreement with data for values of the 
matter density which are bounded from above: Oom < 0.45 (Gold) and f2o m < 0.35 (SNLS). The 
DGP model is excluded at 3a by SNLS and at la by the Gold dataset. Braneworld models with 
future 'quiescent' singularities (at which the Hubble parameter and the matter density remain 
finite but higher derivatives of the expansion factor diverge) are excluded by both datasets. 
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INTRODUCTION 



One of the most remarkable discoveries of the past decade has been the observation 
that the expansion of the universe is speeding up rather than slowing down. The case for an 
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accelerating universe was first made on the basis of high redshift type la supernovae 
and has since received support from deeper and better quality SNe data 0, 0, 0, 
well as observations of the cosmic microwave background and large scale structure 

Theoretically, an accelerating universe can be constructed in a number of distinct ways 
(see [12] and references therein). However three approaches have received considerable at- 
tention in the literature, these are: 



The cosmological constant. The acceleration of the universe is caused by the cosmolog- 
ical constant which satisfies = Ag^ and, hence, has an equation of state p = —p. 
The combination of a A-term and cold dark matter results in the ACDM model which 
appears to provide excellent agreement with cosmological observations and, when com- 
bined with inflationary predictions of an (almost) scale invariant spectrum of density 
perturbations, comprises what may be called the 'standard model of cosmology'. 

However, despite its long and chequered history (since its inception by Einstein in 
1917), a firm theoretical basis for a small A-term has so far eluded researchers 
Indeed, the value of the cosmological constant predicted by quantum field theory is at 
least 10 55 times larger than its observed value p vac = A/8nG ~ 10~ 47 GeV 4 , indicated 
by recent observations. This fact, taken together with the unevolving nature of A, 
suggests that the present epoch may be quite special since Q\ — 2f2o m . The resulting 
cosmic coincidence and the high degree of fine tuning associated with a small A-term 
have lead physicists to look for alternatives to the cosmological constant hypothesis. 

Dark Energy. The expansion of the universe is governed by the field equations of 
general relativity (GR), but one (or more) components of 'matter' violate the strong 
energy condition (SEC) p + 3p > thereby causing the universe to accelerate. To 
this category belong quintessence models, the Chaplygin gas, topological defects and 
numerous other models of dark energy. In general, dark energy can be characterized 
by an equation of state w = p/p, where the observationally determined value of w 
can be used to constrain parameters of a given dark energy model. Although most 
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DE models have w 7^ constant some might say that the cosmological constant with 
w — — 1 also belongs to this category. 

Geometric approaches to acceleration. The late-time acceleration of the universe is 
caused by a departure of space-time physics from standard GR on large scales and/or 
at late times. An important example of this class of models is braneworld cosmology 
according to which our three dimensional universe is a lower dimensional 'brane' em- 
bedded in a higher dimensional 'bulk' space-time. Braneworld models may provide a 
low energy manifestation of string/M-theory lj,ll5|. Within the cosmological context 
braneworld models provide exciting new possibilities some of which are summarized 
below (also see 

(i) 'Quintessential Inflation' [itJ based on the Randall- Sundrum (RS) model may 
provide a compelling explanation of both early and late-time acceleration within a 
single unified setting [18| . (ii) The Dvali-Gabadadze-Porrati (DGP) model can lead to 
an accelerating universe without the presence of a cosmological constant or some other 



form of dark energy [19, |20|. (iii) A family of braneworld models |22j which unify the 



approaches of RS and DGP allow the effective equation of state of dark energy to be 
'quintessence-like' w > — 1 as well as 'phantom-like' w < — 1. In a subclass of these 
models the acceleration of the universe is a transient phenomenon, which gives way to 
matter dominated expansion in the future. The absence of horizons in a transiently 
accelerating space-time is an attractive feature of this scenario since it can reconcile 
current observations of acceleration with the demands of string/M-theory j^j]. (iv) 
Another aspect of the braneworld cosmology 22] is the possibility of fundamentally 



new cosmological behaviour (loitering |24| & mimicry |25() at moderately high red- 
shifts. Loitering and mimicry models remain close to ACDM in the future (hence 
providing excellent agreement with SNe data) while departing from ACDM-like be- 
haviour in the past. The older age of these braneworld models might make them better 
equipped to explain the existence of high redshift QSO's and also allow for a lower 



redshift of reionization than that predicted in ACDM cosmology 
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Whether the increasingly large number of low and high redshift observations can be 
accommodated within the braneworld paradigm is an important subject demanding exten- 
sive exploration. Our purpose in this paper will be more modest, we shall examine the 



braneworld models proposed in (2^ in the light of the Gold SNe data set [7] and the 71 new 
SNe discovered by the Supernova Legaey Survey »J. We shall use this data in ooujuuotiou 
with the recent discovery of the baryon acoustic peak in the Sloan Digital Sky Survey [30( to 
place constraints on the parameter space of accelerating braneworld models. An outline of 
our paper is as follows: in section II we briefly describe braneworld cosmology, while section 
III is devoted to testing braneworld against observations. Our results and conclusions are 
presented in section IV. 



ACCELERATING BRANEWORLD UNIVERSE 



The equations of motion governing the braneworld can be derived from the action [31], |3 



S = M 6 



[ (R 5 - 2A b ) - 2 / k}+ [ (m 2 R 4 - 2a) + / L (h af3 , 0) . (1) 

J bulk ./brane J ./brane ./brane 



Here, R$ is the scalar curvature of the metric g a b in the five-dimensional bulk, and R4 is the 
scalar curvature of the induced metric h a p on the brane. The quantity K = K a ^h a ^ is the 
trace of the extrinsic curvature K a p on the brane defined with respect to its inner normal. 
L(h a p, <fi) is the four-dimensional matter field Lagrangian, M and m denote, respectively, the 
five-dimensional and four-dimensional Planck masses, Ab is the bulk cosmo logical constant, 
and u is the brane tension. Integrations in are performed with respect to the volume 
elements on the bulk and brane. 

The action presents a synthesis of the higher-dimensional ansatzes proposed by Ran- 
dall and Sundrum and Dvali, Gabadadze, and Porrati An important role in 



is played 



paper 



by the m J R4 term. This term was first introduced by Sakharov in a seminal 



2l| to describe the backreaction of quantum fluctuations of matter fields (which, in 
our case, reside on the brane). Its presence is crucial in making the braneworld accelerate 
since it introduces an important length scale 

t = 2m 2 /M 3 . 

On short length scales r £ (early times) one recovers general relativity, whereas on large 
length scales r ^> £ (late times) brane-specific effects begin to play an important role, leading 
to the acceleration of the universe jlH Issj] . 

The cosmological evolution of a spatially flat braneworld is described by the Hubble 



4 



parameter 

4 9 F / 7~A A~! T 7 ^ 

(2) 



where 



and the 'dark radiation' term C/a 4 describes the projection of bulk degrees of freedom 
onto the brane. (Note that the four-dimensional Planck mass m is related to the effective 



Newton's constant on the brane as m = 

The two signs in (J2J) correspond to the two distinct ways in which the brane can be 
embedded in the higher dimensional bulk. Three limiting cases of our model may be of 
interest to the reader: 

1. m = in (J2J corresponds to the well known FRW generalization of the RS scenario 

H 2 , « _ A b ,C (p + a) 2 , . 



2. M = in © gives rise to ACDM 



H 2 (a) = ^ + B. (5) 



3. Finally, by setting a = and Ab = in (J2J) we recover the DGP model [rgj ]. 
The Braneworld models proposed in fall into two main categories: 

• Brane 1 (Bl): The lower sign in (J2J) leads to the following form of the Hubble param- 
eter Q: 

tf 2 (a) = 4 + A eff . (6) 
cr 

The effective cosmological^ 'constant' A e g is composed of two terms: a constant A-term 
and a 'screening term' 



a-=<*+£>4v i+p (? +b -t) 



^ ^ (8) 

A Screening term 
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Since the screening term decreases with time, the value of the effective cosmological 
constant A e g increases. In this respect Brane 1 resembles phantom cosmology which 
has w < — 1. It is important to note, however, that in the braneworld case there is no 
violation of the weak energy condition and also no future 'big rip ' singularity. Indeed, 
from (0) it is quite clear that the universe evolves to ACDM in the future. 

Since our main desire in this paper will be to test braneworld models against ob- 
servations, it will be helpful to recast Eq. (J2J) with the lower sign in the following 
form 

H\z) 

m 



n 0m (i+zy + + 2n £ - 2Jn e Jn 0m (i+z) s + n a + n e + n Ah , (9) 



where z = a /a(t) — 1 is the cosmological redshift, while 

fiom = 3JW' fiCT = 3^br nt= m$' fiAb = _ 6^' (10) 

are dimensionless parameters whose values must be determined from observations. Vt a 
is determined by the constraint relation 



n 0m + n a - 2y/n e y/i + n Ab = 1. (n) 

The difference between ACDM and braneworld cosmology is brought about by the un- 
derlined terms in J^j & Ml\) . 

Brane 2 (B2): The upper sign in (J2J) results in 
H 2 (z) 



H, 



2 



n 0m (i+zy + n a + 2n t + 2Jn t Jtt 0ra (i+zf + n a + n e + n Ab , (12) 



where fie < 1 + fl A and fl a is determined from 



n 0m + n a + 2^n i ^/i + n Ab = 1. (13) 

The difference between ACDM and braneworld cosmology is brought about by the un- 
derlined terms in Alty) & 

The two models Brane 1 and Brane 2 are complementary and reflect the two distinct 
ways in which the brane can be embedded in the bulk. 

From Q & (|12|) it is easy to see that both braneworld models approach the standard 
matter dominated universe at early times [with a small correction term ~ (1 + z) 3 ^ 2 ]. 



At late times the behaviour of the braneworld can differ from both ACDM and SCDM. 
This feature makes braneworld models testable and allows the braneworld scenario to 
provide a new explanation for the observational discovery of an accelerating universe. 

The expansion of the braneworld can be characterized by the deceleration parameter 

q(z) = ^(l+z)-l, (14) 
and the effective equation of state 

2q(z) - 1 

mW = 3[i-n.wr (15) 

From Q, ([12)1 . (|T4"|) and (fi3j) it is easy to obtain the following expression for the 
current value of the effective equation of state 



. ; ^Om / fie , . 

Wo = -l±- — \ - — — - - , (16) 



1 — ^Om V ^Om + Q<t + + 



and we immediately find that Wq < — 1 when we take the lower sign in ([16)1 . which 



corresponds 



41] to Brane 1. The second choice (Brane2) gives wq > —1. 



Mimicry models. It is interesting to note that for values of z and Q\ & &e satisfying 



n 0m (i+^) 3 < Ji + n AbT ^ , (17) 
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Eqs. © and (JT2J) reduce to 



^ CDM (l+^) 3 + l-^ CDM , (18) 



H 2 



where the new density parameter f2^ CDM is defined by the relation 



n ,c D M = a a= Vi±ih. (19) 

a =F 1 v ' ^ 

The braneworld therefore displays a remarkable property called "cosmic mimicry": at 
low redshifts, the Brane 1 universe expands as ACDM (JEJ with f2^ CDM < fi m [Q^ CDM 
is determined by ()19j) with the lower ( "+" ) sign] . The Brane 2 model at low redshifts 
also mimics ACDM but with a larger value of the density parameter f2^ CDM > Q m 
with f2^ CDM being determined by (|TH|) with the upper ("— ") sign. 
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The range of redshifts over which cosmic mimicry occurs is given by < z < 2 m , 
where 

,2/3 

- 1 . (20) 



n l/3 

As (lomonsl rated in 2~| the Hubble parameter in mimicry models departs from that 
in ACDM at intermediate redshifts (z > z m ~ few). This could lead to interesting 
cosmological effects since the age of the high redshift universe in a mimicry model 
can be greater than that in ACDM while the redshift of reionization can be lower 
Since the mimicry models and ACDM are virtually indistinguishable at lower 
redshifts, both are expected to fit the SNe data (at z < 2) equally well. 



COMPARING BRANEWORLD MODELS WITH OBSERVATIONS 



In this paper we shall compare the braneworld model |22j against three sets of observa- 
tions. We briefly summarize each of the data sets which we shall use before proceeding to 
give the results of our comparison. 

1. The Gold SNe data set : As recently as 2003, the entire supernova dataset from the 
two different surveys - Supernova Cosmology Project (SCP) and High z Supernova 
Search Team (HZT), along with low redshift supernovae from Calan-Tololo Supernova 
Search (CTSS) comprised of a meager 92 supernovae , with very few at high 

redshifts, z > 0.7. The method of data reduction for the different teams was also 
somewhat different, so that it was not possible to use the supernovae from the two 
datasets concurrently. The picture changed somewhat dramatically during 2003-2004, 
when a set of papers from both these teams 0, si |f| presented a joint dataset of 
194 supernovae which used the same data reduction method. This new data resulted 
in doubling the dataset at z > 0.7. Not all these supernovae could be identified 
beyond doubt as Type la supernovae however, in many cases complete spectral data 
was not available. In early 2004, Riess et a/.j7j reanalyzed the data with somewhat 
more rigorous standards, excluding several supernovae for uncertain classification or 
inaccurate colour measurements. They also added 14 new high redshift supernovae 
observed by the Hubble Space Telescope (HST) to this sample. The resultant sample 
comprises of 157 supernovae (the furthest being at redshift z = 1.75) which have been 
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classified as Type la supernovae beyond doubt- the 'Gold' dataset. We shall be using 
this 'Gold' dataset as our first supernova sample. 

2. The Supernova Legacy Survey SNe data set (SNLS) : The SuperNova Legacy Survey ^] 
is an ongoing 5-year project which is expected to yield more than 700 spectroscopica 
confirmed supernovae below redshift of one. The first year results from this survey 
have provided us with 71 new supernovae below z — 1. We shall use these 71 SNe 
together with the already available low-z supernova data, ie a total of 115 SNe, as our 
second supernova sample. 

3. The Baryon Acoustic Oscillation Peak (BAO) : A remarkable confirmation of the 
standard big bang cosmology has been the recent detection of a peak in the correlation 
function of luminous red galaxies in the Sloan Digital Sky Survey (3^. This peak, 
which is predicted to arise precisely at the measured scale of 100 h _1 Mpc due to 
acoustic oscillations in the primordial baryon-photon plasma prior to recombination, 
can provide a 'standard ruler' with which to test dark energy models. Specifically, we 
shall use the value [^l 

A v 7 ^ r 1 dr. 12/3 



■ 1 


r z i dz ' 


- Zi 


Jo h{z) . 



0.469 ±0.017, (21) 

h(zi) l ' 6 L z\ Jo h{z) . 
where h(z) = H(z)/Hq is defined in (Q and (|12p for Branel and Brane2 respectively, 
and Z\ = 0.35 is the redshift at which the acoustic scale has been measured. 



Methodology and Results 

For the supernova data, we shall use the \ 2 minimization where 

^(u a n \ sr \yfoA z i> Hq> ^om» Pj) - yi? / 99 x 

i °i 

Here, y^ is the data at redshift of z^ and (7j is the uncertainty in the individual yi, and pj 
are the braneworld parameters (fl^,fl\ b for Branel and Brane2). For the 'Gold' dataset, 
Vi = A*o,i = m B — M = 5\ogdL + 25, the extinction corrected distance modulus for SNe at 
redshift Z{. The error er(/io,i) includes the uncertainty in galaxy redshifts due to a peculiar 
velocity of 400 km/s. For SNLS, yi = fiB,i = itib — M + a(s — 1) — (3c = 5 log 10 (<iz / /10pc). 
The error o~(fiB,i) includes effects due to a peculiar velocity of 300 km/s. We assume a flat 
universe for our analysis. 
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We note that, for the SNLS data, we have to deal with two additional parameters a, (3 
during the minimization. However, the error a is dependent on a, (3. Therefore, if we 
minimize with respect to all the parameters, the process will be biased towards increasing 
the errors in order to decrease the x 2 ■ To avoid this, we refrain from varying a,f3 together 
with the cosmological parameters. At each iteration a, (3 are fixed, while the remaining 
parameters are varied to obtain the minimum x 2 , then the values of a, (3 are changed for the 
next iteration. This process is continued till the global minimum is obtained. This method is 
equivalent to the method followed by the SNLS team, and for flat ACDM, our results concur 
with those in gj]. In presenting the final results, the nuisance parameters a, (3, M (where M 
depends on H ) are marginalized over to obtain bounds on the cosmological parameters of 
interest. 

For Brane 1 and Brane 2, the cosmological parameters to be estimated are 
n 0m , Qi and Q\ b (fl a is calculated from Eqs (fTTj) and ()13|) respectively for Brane 1 & Brane 2). 
After marginalizing over the statistical nuisance parameters, we obtain the three-dimensional 
probability distribution in the (f2om, ^A b ) space: P(f2 0m , ^A b )- We perform maximum 
likelihood analysis on the system with the priors < Oom < 1, &e > 0, Q\ b > 0. For 
Brane 2, we use the added constraint < 1 + ^A b . For Brane 1 & Brane 2 the constraint 
relations ()13|) combined with Q K = 0, set the lower bound Q\ b > — 1. However since 
^A b > is a more physically appealing model (it includes anti-de Sitter space (AdS) bulk 
geometry), we choose this as a prior for further analysis. 

We may add further information to the analysis from the baryon acoustic oscillation 
(BAO) data. We obtain the joint probability distribution for the SNe data and the BAO 
data as -P(Xsn c +bao) = P{Xstie)P{A), where A is the quantity defined in eq (j3J) and we 
assume that it follows a Gaussian probability distribution with mean A = 0.469 and an 
error of o = 0.017. 

In the figure ^ we show the results for Brane 1 and Brane 2 using both Gold and SNLS 
data, in conjunction with the baryon acoustic oscillation peak (BAO). We show the reduced 
X 2 per degrees of freedom as a function of flom, marginalized over f^, flA b , Hq. For Brane 1, 
we find that the supernova data alone, in both cases, favours a somewhat larger value of O 0m 
at the minimum, with 'Gold' preferring a higher value than SNLS. When used in conjunction 
with BAO, however, both datasets prefer a matter density of f2om ~ 0.26. For Brane 2, the 
preferred value of O 0m is around f2 0m = 0.2 for SNLS and around f2 0m = 0.3 for 'Gold'. 
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TABLE I: Best-fit ilo m — f2g with corresponding la errors for the Branel model, with the present value of 
the equation of state Wq and the best-fit \ 2 for the two supernova datasets. The best-fit \ 2 for ACDM is 
also shown for comparison. 
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TABLE II: Best-fit f^om — &i with corresponding la errors for the Brane2 model, with the present value 
of the equation of state wq and the best-fit \ 2 f° r the two supernova datasets. The best-fit x 2 f° r ACDM is 
also shown for comparison. 
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When used with BAO, once again, fiom — 0.26 is preferred. 

During the analysis we find that for both Brane 1 and Brane 2, the results are very weakly 
dependent on the bulk cosmological constant ^A b , and marginalizing over Q^ b does not affect 
the results very much. In further analysis, we therefore fix the value of Q\ b to its best-fit 
value of ^A b = 0. 

In the table|l]we show the best-fit Qom — &e with la errors around it for the Brane 1 model 
for the two supernova datasets, with the x 2 a t the best-fit. The best-fit \ 2 for ACDM is also 
shown for comparison. We know that the supernova data alone is not able to place strong 
constraints on the value of the matter density Q 0m for Brane 1 models for both datasets 
(see upper panel of figure EJ). If we use the baryon acoustic oscillation peak in conjunction 
with the supernova data, then there are stringent and realistic constraints on the value of 
ilom, therefore we show the x 2 f° r t ne joint probability distribution only. For the Brane 1 
model, Qg cannot be negative, and wq < — 1 always, and we show the bounds on and 
Wq accordingly. For the Gold dataset, the supernova and BAO data in conjunction favour 
a ACDM universe with Q 0m = 0.2741°,;°,! and Q t = 0.0+ 031 . The "effective" equation of 
state at present is w = — 1-0_ 0057 . For the SNLS dataset, a braneworld model with a small 
value of Q e is slightly favoured over ACDM, with Q 0m = 0.2691^1 and Q e = 0.05418^. 
The present value of the "effective" equation of state is wo = — 1.073i°)]o35- Since Vt^ = 
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FIG. 1: The reduced x 2 P er degree of freedom as a function of Oo m , marginalized over f^, £l\ b , Hq 
for Brane 1 (top) and Brane 2 (bottom). The left panel shows results for the Gold Supernova data 
while the right panel shows results for the SNLS data. The dashed line is each panel is obtained 
by using supernova data alone, while the solid line uses both SNe data and the baryon oscillation 
peak. 
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FIG. 2: The redshift variation of the "effective" equation of state w(z) for the Brane 1 (top panel) 
and Brane 2 (bottom panel) models using the Gold+BAO (left panel) and SNLS+BAO (right 
panel) datasets. The light grey contours denote the la errors around the best-fit, and the dashed 
line represents ACDM, which is the upper (lower) limit for w(z) for the Brane 1 (Brane 2) model. 
In the top right panel, the thick solid line represents the best-fit for the SNLS+BAO data for 
the Brane 1 model. For the other three cases, the best-fit is at the ACDM line. We see that the 
behaviour of the "effective" equation of state of the braneworld models can be markedly different 
from ACDM within la even for the small values of Qe allowed by the data, especially for Brane 1. 
For the Gold data, which extends to higher redshifts, we can even see the existence of a pole in 
the "effective" equation of state for the Brane 1 models. 

13 



represents ACDM, this implies that BAO and SNe data together choose a best-fit Brane 1 
model that is either equivalent to, or very close to, the ACDM model. However, within 
la, the value of Qg allowed is sufficiently large to allow for interesting braneworld behavior, 
especially for the SNLS data. 

In the table |H1 we similarly show the best-fit f2 0m — fie with la errors around it for 
the Brane 2 model for the two supernova datasets, with the x 2 & t the best-fit, along with 
the best-fit for ACDM. As before, since for the Brane 2 model, fie cannot be negative, and 
Wo > — 1 always, we show the bounds on Qi and Wq accordingly. In this case, the current data 
appears to favour a ACDM universe over the Brane 2 model for both datasets. Somewhat 
larger values of fie are allowed for Brane 2 within la for the Gold data than for the SNLS 
data. 

The low values of fie allowed by the current data for both the braneworld models can 
still give rise to sufficiently interesting behaviour. To demonstrate this, in figure 121 we show 
the la confidence levels for the "effective" equation of state w for the Brane 1 and Brane 2 
models using Gold+BAO and SNLS+BAO data. We see that, even for the small values of fie 
allowed by the data, the "effective" equation of state looks quite different from the w = —1 
cosmological constant model within la especially for the Brane 1 model. In fact, for the 
Gold data, which extends to higher redshifts, we can see the existence of a singularity in the 
"effective" equation of state at z ~ 1.6 for the Brane 1 model. We would like to emphasize 
here that the presence of such poles in the "effective" quantity w(z) does not signal to any 
inherent pathologies of the braneworld models described here, since the scale factor and its 
derivatives remain well behaved throughout the evolution of the universe. Indeed, the reason 
for the occurance of a pole is simple and has to do with the fact that the density parameter 
fl m (z), which increases with increasing redshift, crosses unity at high z 22J,|24|. This results 
in a pole in the effective equation of state, since w(z) in ()15|) diverges when fl m (z) — > 1. It 
should be stressed that the diverging equation of state is a signature of these braneworld 
models. Although it is difficult to reconstruct a diverging equation of state using standard 
parametrizations of dark energy, the equation of state of such braneworld models can be 
successfully reconstructed using other reconstruction techniques such as smoothing of the 
supernova data |34j. So, it is hoped that with better quality data, it should be possible 
to test braneworld models by studying the behaviour of the effective equation of state at 
high z. At present we see that, within la, braneworld models that are distinctly different 
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from the cosmological constant are able to satisfy the current supernova data and therefore 
remain a possible candidate for dark energy along with the cosmological constant. 

In figure |2J we further explore the f^om — &e parameter space for the Brane 1 and Brane 2 
models. We see from the la, 2a, 3a contours that Brane 1 satisfies the SNLS dataset for 
a larger region in the f2 0m — parameter space as compared to the Gold dataset. BAO 
does not depend very strongly on the value of Qe, rather, it is sensitive to the value of f2 0m - 
Therefore, in conjunction with BAO, Brane 1 provides a better fit to the SNLS data than 
to the Gold data. The situation is just the opposite in the case of Brane 2, which provides 
a better fit to the Gold data as compared to the SNLS data. 

An interesting feature of braneworld models is that for a finite region of parameter space, 
the Brane 2 universe can expand towards a 'quiescent' future singularity at which the en- 
ergy density and the Hubble parameter remain well behaved, but higher derivatives of the 
expansion factor (a, a etc.) diverge j^j]. From (JT2jl we see that a braneworld model which 
satisfies 

n a + n t + n Ab < o. (23) 

will run into a future singularity at the redshift 

*,= (- ii ° + "' + "M ' /3 -l. (24) 

The time of occurrence of the singularity (measured from the present moment) can easily 
be determined from 

dz 



r . = t( , = *,_ t( , = )=/, rT _, (25) 



where H(z) is given by (|12|). 

From the lower panel of figure El we see that universes which terminates in a 'quiescent' 
future singularity are excluded at the 3a confidence level for both SNLS and Gold datasets 
when used in conjunction with the baryon oscillations. 

As noted in section „ the DGP model forms a special case of our braneworld cosmology 
obtained by putting = Q a = in equation H21 We see that using the SNLS data together 
with BAO, we may narrowly rule out the DGP model of braneworld dark energy at 3a (thick 
solid line in right lower panel of figure EJ). However, for the Gold data and BAO, the flat 
DGP model is acceptable within 2a (thick solid line in left lower panel of figure l3|) [d^ . 
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FIG. 3: Confidence levels in the f^om — &i plane for Brane 1 (top) and Brane 2 (bottom). Hq is 
marginalized over and 0,\ h is fixed at the best fit value of £l\ b = 0. The left panel shows results 
for the Gold Supernova data while the right panel shows results for the SNLS data. The solid 
lines in each panel represent the 1, 2, 3a confidence levels obtained by using supernova data alone, 
while the dashed lines are the 1,2,3a contours for the baryon oscillation peak. The light grey, 
medium grey and dark grey contours represent the 1,2,3a confidence levels when the supernova 
data is used in conjunction with the baryon oscillation peak. The region to the right of the thick 
dotted line in the lower panels represents Brane 2 models which will undergo a future 'quiescent' 
singularity. The thick solid line in the lower panels represents the DGP model. 
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We therefore conclude that while the flat Brane 1 and Brane 2 models are able to satisfy 
the Gold and SNLS data data over a reasonable region of parameter space, both the DGP 
model and the model with a quiescent future singularity are in tension with the data. For 
both the supernova datasets, the more general flat Brane 1 and Brane 2 models are able to 
satisfy the data over a reasonable region of parameter space. 



CONCLUSIONS 

We have demonstrated that high redshift type la supernova data when combined 

with the recent discovery of baryon oscillations in the SDSS [3(| can serve to place significant 
constraints on the parameter space of braneworld models. 



Our results for the Gold data set are in broad agreement with the earlier work of 
who tested braneworld models against an early SNe data set (also see j^Zl). Our results for 
the Supernova Legacy Survey (SNLS) are in good agreement with those of J3| who recently 
tested the DGP model using SNLS and baryon oscillations. Since the DGP model forms a 
subset of the braneworld models analyzed by us we find that the SNLS data together with 
baryon oscillations rule out this model at 3 a. However we also find that the DGP model 
is more strongly constrained by SNLS than by the Gold data set of j7j, which allows the 
DGP model at 2a. Our analysis also shows that baryon oscillations in conjunction with SNe 
data rule out a class of braneworld models in which the universe encounters a 'quiescent' 
future singularity where the density, pressure and Hubble parameter remain finite but higher 
derivatives of the scale factor diverge. 

Our analysis indicates that the Gold and SNLS supernovae place slightly different con- 
straints on the braneworld parameters. Thus although figures ^ 01 clearly show that the 
Braneworld models analyzed by us agree well with both sets of SNe data, the Gold data set 
accommodates larger values of Qom > 0.25 than SNLS (ffom > 0.2) for the Brane 1 model. 
In the case of Brane 2 smaller values f2o m <, 0.35 appear to be favoured by SNLS than by 
Gold (fi m <, 0.45). Thus Brane 1 models fit better to the SNLS data, while Brane 2 models 
fit better to the Gold dataset. 

In this connection it is interestin g to note that the recent analysis of evolving dark energy 
models using Gold and SNLS data [39j found somewhat different trends in the evolution of 
dark energy favoured by these two data sets. It is hoped that improvements in the quality 
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and quantity of future SNe data will allow tighter constraints to be placed on dark energy 
models. We would like to draw attention to the fact that in this paper we have explored 
the region of parameter space in which is reasonably small. The reason for this is that 
performing a likelihood analysis for the full (formally infinite) region of parameter space is 
computationally very expensive and so we have restricted our analysis to a finite region of 
parameter space corresponding to ^A b 1, ,$ 1 which includes ACDM and the DGP 
braneworld as subclasses. In this context it is interesting to note that braneworld models 
with fundamentally new properties can arise for large values of Q\ b ^ 1- For instance, 
the mimicry model briefly touched upon in section II, mimics ACDM at low redshifts while 
departing from the latter at higher redshifts. The presence of a 'dark radiation' term (which 
quantifies the projection of higher dimensional 'bulk' effects on to the brane) allows for even 
more radical departures from standard cosmology by permitting the universe to 'loiter' at 



moderately large redshifts 1 <C z <C 1000. Both 'loitering' |24| and 'mimicry' 25] models 
predict a lower redshift of reionization and a longer age for QSO's and other high redshift 
objects when compared with ACDM. However at low z both models remain very close to 
ACDM and for this reason are unlikely to be distinguished from the latter on the basis of 
SNe data alone. For this reason we have not included 'loitering' and 'mimicry' models into 
the present analysis but hope to return to these models in a companion paper. 
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